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The analytical theory for the voltage-current characteristics of the large inductance (L > lOOpH) 
high- Tc DC SQUIDs that has been developed previously is consistently compared with the computer 
£S| . simulations and the experiment. The theoretical voltage modulation for symmetric junctions is 

shown to be in a good agreement with the results of known computer simulations. It is shown 
that the asymmetry of the junctions results in the increase of the voltage modulation if the critical 
current is in excess of some threshold value (about 8/j.A). Below this value the asymmetry leads to 
^N^j ■ the reduction of the voltage modulation as compared to the symmetric case. The comparison with 

the experiment shows that the asymmetry can explain a large portion of experimental values of the 
voltage modulation which lie above the theoretical curve for symmetric DC SQUID. It also explains 
experimental points which lie below the curve at small critical currents. However, a significant 
portion of these values which lie below the curve cannot be explained by the junction asymmetry. 
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I. INTRODUCTION 



The new type of superconductors, which has been discovered in the end of the last century by Bednorz and Muller, 
is widely used in the modern SQUID systems. The majority of high Tc DC SQUIDs are based on the YBCO thin 
films and have the different types of designiSi^. However, the adequate theory of the voltage-current characteristic 
(VCC) of the high- Tc DC SQUID, which would predict its transfer function and energy resolution, still not exists. In 
recent time intensive computer simulations and theoretical studies have been performed to investigate the dependence 
• of high- Tc DC SQUID behavior on various factors^£*&, but a marked disagreement of the numerical simulations 
with experiment is still observed: the experimental transfer functions in many cases are much lower than the values 
predicted by theory and computer simulations; the white noise is about ten times higher than predicted. This is 
one of the most important unsolved problems, which seriously hinders the optimization of high Tc DC SQUIDs for 
i applications. 

^C) ' Up to now the high-Tp DC SQUIDs have the significant parameter dispersion, but the reasons of such dispersion 
are not established. One of the possible reasons for the dispersion could be attributed to the junction asymmetry of 
SQUID interferometer (unequal critical currents or (and) normal resistances), which for grain boundary junctions is 
| about 20-30 percents due to on chip technological heterogeneity. 

Recently, the theory of the voltage-current characteristic of the high Tc DC SQUID, which expands the validity 
range of the Chesca's analytic theorj£ on the DC SQUIDs with large inductance L > 100 pH and accounts both 
for the symmetric and asymmetric DC SQUIDs, was developed^. The theory is based on the perturbation solution 
of the two-dimensional Fokker-Planck equation (2D FPE) that describes the stochastic dynamics of DC SQUID in 
I . presence of the large thermal fluctuations. 

In the paper the results of the analytic theory of VCC of high Tc DC SQUIDS are compared with the computer 
simulations and with the experiment. It is shown that the theoretical voltage modulation for symmetric junctions 
is in a good agreement with the results of known computer simulations. It is also shown that the asymmetry of 
the junctions results in the increase of the voltage modulation if the critical current is in excess of some threshold 
value (about 8/iA). Below this value the asymmetry leads to the reduction of the voltage modulation as compared 
to the symmetric case. The comparison with the experiment shows that the asymmetry can explain a large portion 
of experimental values of the voltage modulation which lie above the theoretical curve for symmetric DC SQUID. It 
also explains experimental points which lie below the curve at small bias currents (less about than 10/xA). However, a 
significant portion of these values which lie below the curve at larger bias currents cannot be explained by the junction 
asymmetry. 

The paper is organized as follows. In the section II we present in a concise form the main theoretical results for 
symmetric and asymmetric DC SQUIDs. In Section III we compare theoretical voltage-current characteristics for 
symmetric DC SQUID with the computer simulations of stochastic dynamical equations of DC SQUIDs which have 
been made earlier by other authors. In this section we also study in detail the influence of the junction asymmetry 
on the voltage modulation as compared with symmetric case. In Section IV we compare the theory with experiment 
and show that the junction asymmetry can explain a large portion of experimental points which lie well above the 
theoretical voltage modulation curve for symmetric DC SQUID. 
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II. THE MAIN RESULTS OF THE ANALYTIC THEORY OF VCC OF HIGH T c DC SQUID 

A. Symmetric DC SQUID 

We consider a symmetric DC SQUID with equal critical currents of junctions, Iqx — Ic2 = Ic, equal normal 
resistance, R\ — R2 = R, and loop inductance L. The equations, describing such SQUID, have the following form: 
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where $0 = h/2e is a quantum of magnetic flux, $ is a magnetic flux trapped in the interferometer loop, &x is 
external magnetic flux, I is a bias current, <p = 7r<I>/<f>o, 8 — 7rA/<I>o. The quantities 7jv±(t) are independent stochastic 
variables related to the Nyquist current noise of junctions: 



(l N± (t)I N± (t')) = ^f6(t-t') 



(3) 



where k is the Boltzmann constant, T is the absolute temperature. 

The output voltage across SQUID is the low frequency part of the equation @, which is averaged over the noise: 
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Thus, the stochastic equations 0]), are described DC SQUID behavior in presence of the thermal fluctuations. 
The solution of these equations depends on the following parameters: screening parameter (3 = 2LIc/&o, noise 
parameter T = 2TrkgT/^>Qlc and dimensionless inductance a — L/Lp, where Lp — (&q/2it) 2 /kgT is a fluctuation 
inductance, equaled 100 pH at = 77 K. But only two of them are independent due to relation a = 7r/3r. 

The known theoretical approaches to the solution of equations QJ, J2J are based on the analysis of the equivalent 
two dimensional Fokker-Planck equation for the distribution function P($, A): 
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where W is a two dimensional potential energy for DC SQUID: 
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Ej = QqIc/2it is a Josephson coupling energy of two junctions in parallel, i = I /Ic, 4>x = 7r<I>x /$o is a dimensionless 
external magnetic flux. 

At the present time there exist some analytical solutions of equation (J5J, which are valid at different ranges of 
0,T,a. 

I) The exact analytical solution of FPE can be obtained in the small inductance limit (L — > 0, j3 « I, a << I). 
In this case SQUID is equivalent to a single Josephson junction with normal resistance R/2 and critical current 
27(7 1 cos 4>x I ■ The voltage across such single junction is well known and can be obtained from Ambegaokar-Halperin 
formS: 
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W(x) = (i/T) x + (2/T) cos ^ x cos x (9) 

2) The approximate analytical solution of FPE JSJ was obtained ir£. This solution is valid at [3 < 0.3 and relatively 
high thermal fluctuations level (r > 1). 

3) The original method for the solution of equation JS} for SQUID with large inductance (a > 1) was presented 
The method is based on the perturbation expansion of the solution of 2D FPE (JSJ over small parameter 
e = cxp(— a/2). The result is the following expressions for voltage V across a SQUID and voltage modulation 
AU = V{<t> x = tt/2) - V(0 X = 0): 

1 " 7fl iexp(-a/2)cos(2^)/(i,r) (10) 



RI C p(i, 2r,0) 2 



AV 

— =exp(-a/2)/(i,r) (11) 



where the value 2r, 0) was defined in (|SJ|. 
According to& 



where 
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I n is a modified Bessel function. 

These expressions are valid for a > 1 and any values of /3, and T which are consistent with the condition a — Trf3T. 
Therefore, they can be applied for the analysis of a majority of practical high Tc DC SQUIDs with T ss 0.05 — 1, 
ft > 1, ot > 1. However, it should be remembered that Eqs.Q and (|llfl are the approximate expressions which 
account for the first order term in the perturbation expansion of the voltage over small parameter e = exp(— a/2). 



B. Asymmetric DC SQUID 

We describe the junction asymmetry in terms of the asymmetry parameters 7 and p, which are defined according 
to: Ici = {l + l)Ic, Ic2 = (l - i)I c ,Ri = R/0- + P),R2 = ity(l-p), where 

Ici+Ic2 Id - IC2 ,n 

/C = ^2^ ; ^TcTTfe (14) 

R — ^^^ 2 p — ^ 2 ^ > (15) 
Ri + R2 Ri + R2 

The perturbation method, developed in£, has been applied to obtain the expressions for the voltage and its modulation 
across an asymmetric SQUID with large inductance^. Corresponding expressions have the following forms: 

1 li ' 87r 3 e- a / 2 — {Scos(2(px + t )+Qsm(2<fx+t )} (16) 
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The quantities S 1 and Q have the forms: 
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For a > 1 the expressions (|16|I - H24|I are valid at any values of asymmetry parameters p and 7 and at any values of (3 
and r, which are consistent with the condition a = 7r/3r. 

III. VOLTAGE-CURRENT CHARACTERISTICS 
A. Symmetric DC SQUID 

According to l|10(l . the influence of SQUID inductance on VCC appeared as the reduction of the apparent value of 
the critical current (Fig. . Such behavior is due to the the suppression and masking of the critical current by the 
significant thermal current fluctuations in the interferometer loop. The more the inductance the more the suppression 
of the critical current. If a >> 1 the second term on the righthand side of the expression (|10fl may be neglected and 
from comparison of i|l(J|) with (0J we see that in this limit the DC SQUID is equivalent to a single Josephson junction 
whose critical current is twice as less as that for the case L — 0. 

There are two DC SQUID parameters which are easily measured: the bias current, /, and the voltage modulation, 
AV. By the tuning of the bias current the value / = Imax, which corresponds to the maximum of the voltage 
modulation AVmax, can be found. On the practice this guarantees the maximum of the SQUID transfer function 
dV/d^x- 

As is seen from (|11[) its righthand side is the product of two terms , one of them exp(— a/2) depends on the 
SQUID inductance only, the other one depends on the bias / and critical Ic currents. The first factor describes the 
suppression of the critical current by the noise current in the interferometer loop. This is similar to the suppression 
effect in the interferometer loop with a single Josephson junction 15 . The second factor describes the critical current 
suppression by the thermal fluctuations, which is similar to the suppression effect in a single Josephson junction 9 . 
The factorization allows us to carry the first factor to the left hand side of expression (|ll(l , so that we consider below 
the reduced modulation AVr = exp(a/2)AU/i?, which depends on the critical and bias currents only. 

The typical dependence of the reduced voltage modulation on the bias current I at different critical currents is 
shown in (Fig. [2J. The curves, corresponding to the different Ic's, are shifted along the current axes and for given Ic 
each dependence AVr(I) has a well defined maximum AVr(Imax) = AVr^max a t the corresponding value of the 
bias current. 

From the equation <|ll|l for a set of fixed values of Ic we have computed a set of the values of maximum voltage 
modulation AVr.max with the corresponding values of bias current Imax- In this way we have obtained the table 
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of values AVr^max{Ic,Imax)- With the aid of the table we draw the dependence AVr } max{Ic), which is shown 
in Fig. [3J It is obvious, that this curve gives the upper bound of AV for any symmetric DC SQUID. The different 
points on the curve correspond to the different bias currents I max, at which AV reaches its maximum. 

However, as was mention above, the critical current of high-Tp interferometer cannot be measured directly with 
a good accuracy because of large thermal fluctuations. Therefore, it is useful in practice to use the dependence of 
the maximum modulation signal AVmax on the corresponding value of the bias current Im ax- Such dependence, 
obtained from the table of values AVr,max{Ic, I max), is shown on Fig. 0] The different points on this curve belong 
to the different values of Ic- Every point on the curve is the maximum point on the corresponding curve from Fig. [21 

The characteristic feature of these two curves AV^MAxilc) and AVr^max(Imax) is the saturation of AVr.max at 
large critical (Ic > 50fiA) and bias currents. From the curve shown on Fig. [3|the critical current can be obtained from 
the measured value of AVmax- I n addition, this dependence allows one to predict the maximum voltage modulation 
AVmax, if the critical current of a SQUID is known from the direct measurements. 

The influence of inductance on AVmax at different values of the noise parameter, T is shown on Fig. [5] The 
curves were calculated from Eq. (fTT|l . It can be seen, that the increase of SQUID inductance leads to the reduction 
of the voltage modulation in accordance with the scaling law exp(— L/2Lf)- In addition, the increase of the noise 
parameter, T also leads to the decrease of the voltage modulation. 



B. Critical current 



As is known, it is difficult to measure the critical current of high-Tc Josephson junctions with a good accuracy 
because VCC is washed out by the thermal fluctuations. This problem can be solved, if we relate the critical current 
with bias current Imax, which can be measured directly in SQUID scheme. The idea was realized in the paper—, 
where, based on the numerical simulations of the paper—, the following approximate expression for Imax has been 
suggested: 



Imax = 2I C (l - VW^) (25) 
From Eq. H'25f) the critical current can be expressed in terms of the well measured bias current, Imax- 



which for T = 77K becomes: 



Imax 



0.514 (l + v/f + 1.945/max) (27) 



where Ic, Imax are in /iA 

The Eq. (|25[1 has been obtained for low- Tc junctions^ and applied to high Tc SQUIDs in— Therefore, it is 
interesting to compare (|26[1 with our theory. As is seen from Fig. [3 the expression (|llf> establishes a single- valued 
dependence between the critical current and the bias current. The bias current is the well- measured parameter in 
SQUID scheme. If Imax is known from the experiment the expression (|f I \ permits to find the critical current at 
which the voltage modulation reaches its maximum. From lllfl for a set of critical current values, Ic we have found 
a set of the bias currents Imax, which give the maximum voltage modulation. The comparison of Drung's expression 
(|27H with our theory is shown on Fig. EI As is seen from the figure Drung's expression l|27|) (solid line) and the theory 
(crosses) give approximately similar results. The deviation between Drung's curve and theoretical points is no more 
than 10% in the whole range of Imax's- 

Below we compare the analytical expression (|llfl with the results of the computer simulations of the stochastic DC 
SQUID equations Q and @, obtained by other authors. The normalized voltage modulation AV/RIc as a function 
of the bias current at the different SQUID inductances is shown on Fig. The curves marked by black stars are the 
theoretical ones, calculated from Eq. (|ll|l (the curves for L = 94pT, L = 157pT) and from the Ambegaokar-Halperin 
expressions J7J), © (the curve for L = 0). The curves marked by open circles obtained by computer simulations of 
exact stochastic equations Q and J2J 17 - All calculations and simulations were made for T = 77 K and noise parameter 
r = 1 [Ic = 3.23^, A). As is seen from Fig. [7| the simulated curves are close to the calculated ones. Therefore, our 
expression (fTTI) is a good approximation of the exact solution of stochastic equations and J2J) ■ 
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In addition, we compare our theory with the two well known expressions for the maximum transfer function, 
V$ = dV/ d<& which have been obtained by the empirical fit to the computed values obtained from the simulations of 
the exact stochastic equations of DC SQUID. The first expression is obtained inS and is valid for /3 > 0.5, T < 1: 

AV MAX _ 7.3/30- 



X c R tt(1+/3) (80r/3) ' 4 + 0.35 (4r/?) 2 - 5 

The second one is the widely used expression of EnpukvA 

AVmax 4 ( 2 k B TL\ 

= — ; exp — 3.57T ^— (29 

I C R tt(1 + /3) l \ $g ) y ' 

Here the transfer function was recalculated to AV (AV = V^q/tt) assuming a sine shape of the voltage-to-flux curve. 
The comparison of our theory with expressions (|28|) and (|29[1 for several values of the inductances is shown on Fig.lJS}. 

It is seen, that the theoretical curve and the curve of Enpuku (|29|l reach the saturation at approximately 
Ic > 40/1^4 while the curve of Kleiner l|28|) has a constant non vanishing slope. This slope is probably due to the fact 
that the right hand side of (|28|l is equal, in fact, to the transfer function V$ which has actually been calculated ir£. 
For relatively high critical currents (small T's) the sine approximation (AV = V$$o/tt) we made for the shape of the 
signal is not very good due to the distortion of the signal shape. It is also worth noting that the SQUID inductance 
affects AVmax in different ways. The Eq. (|28|) always gives the highest values, except for a = 2 for Ic < 80/j.A 
(Fig. The Eq. (|29|l always gives the lowest values, except for a = 1 (Fig. For large inductances a > 2 the 
theoretical curve is always higher. The influence of inductance is more pronounced for the voltage modulation given 
by Eqs. (|28|l and (|29p. For example, for Ic = 100/1^4 from a = 1 to a = 2.5 the theoretical voltage modulation is 
reduced by a factor of two, while the reduction factor for the Kleiner's and Enpuku's expressions is six and nine, 
respectively. 



C. Asymmetric DC SQUID 

The numerical simulations of stochastic differential equations, which govern the dynamics of DC SQUID are very 
time consuming for practical high Tc DC SQUIDs due to the large thermal fluctuations and large loop inductance. 
This is why in most cases the investigations of asymmetric DC SQUID are restricted to the computer simulations in 
small inductance limit (a « i^iiiiiiSiiii^. Since our theory is valid for asymmetric SQUIDs with large inductance 
we cannot compare it with the results of other authors obtained for asymmetric SQUIDs with small inductance. 

The practical importance of asymmetric DC SQUIDs is that they have a higher transfer function (the slope of the 
voltage-to-flux curve) as compared to the symmetric casei2*ii*i2iiii4. This property is generally attributed to the 
distortion of the voltage-to-flux curve. The shape of this curve significantly differed from the sine shape that results 
in the high steep of the slope. For large inductance SQUID the voltage-to-flux curve has a sine shape as evident from 
Eq. (|16|l . Therefore, we carefully studied the influence of the junction asymmetry of large inductance interferometers 
on the voltage modulation. We have found that in the large range of asymmetry parameters the voltage modulation 
of asymmetric DC SQUID with large inductance can be significantly higher than that for symmetric SQUIDs. 

In principle, the asymmetry parameters 7 and p are independent of each other and, therefore, they may present 
in the junctions in any combination. From our calculations we have chosen for presentation here only three types of 
the junction asymmetry which generally give a correct picture of how the voltage modulation is influenced by any 
type of asymmetry. Below we present the numerical results for the current asymmetry (p = 0,7 ^ 0), the resistance 
asymmetry (p 7^ 0,7 = 0) and the geometric asymmetry (p = 7 ^ 0) 11 . The effect of different types of the junction 
asymmetry on the dependence of the voltage modulation on the critical current is shown on Fig.|5|and Fig. 1101 As is 
seen from the figures, the general trends of the asymmetry are the increase of the voltage modulation and non vanishing 
slope of the curves as compared with the symmetric junctions. The slope being increased with the increase of the 
asymmetry. For the geometric asymmetry the increase of the voltage modulation is seen in the whole range of the 
critical currents (Fig.^J, while for the current or resistance asymmetry there exists the range of relatively small critical 
currents (approximately less than 10fiA) where the voltage modulation of asymmetric DC SQUID is lower than that 
of the symmetric one (Fig. llOfl . In order to clarify the picture we made a careful study of the dependence of AVr.max 
on the asymmetry parameters at the given values of the critical current. It appears there exists some threshold value 
of the critical current, approximately in the vicinity of 8(iA, which divides the whole range of the critical currents 
in two parts. Below the threshold the current and resistance asymmetry always leads to the decrease of the voltage 
modulation as compared to symmetric DC SQUID fFig. II lb. Fig. II lb). The geometric asymmetry does not change 
the voltage modulation up till approximately p = 7 = 0.5 with the subsequent decrease of the voltage modulation 
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(Fig. ITTb). Above the threshold the curves AVr,max(7j P — 0,7c = const) and AVr,max(7 = ^,p,Ic — const) 
have a clear maximum approximately at the vicinity 7 = 0.25 and p — 0.25, respectively (Tig. 112b. Fig. 112b). This 
maximum lies approximately 1.5 times higher a symmetric value for the voltage modulation for current asymmetry 
and 1.3 times higher a symmetric value for the voltage modulation for resistance asymmetry. It is worth noting that 
in the wide range of the current asymmetry (0 < 7 < 0.7) or of the resistance asymmetry (0 < p < 0.7) the voltage 
modulation is higher than its symmetric value. However, the geometric asymmetry leads to the significant increase 
of the voltage modulation in the whole range of the critical currents (Fig. 112b). 

IV. COMPARISON WITH EXPERIMENT 

In order to compare our theory with experiment we took two groups of DC SQUIDs. The first group comprised 
about 50 SQUIDs which have been chosen before for the same purposei&jiS. All of these SQUIDs are single layer 
ones, using 100 or 200 run thick Y Ba-iCu^Oi^x films deposited by laser ablation onto SrTiO% bicrystal substrates 
with 24° or 30° misorientation angles, both having symmetrical configuration. The technology is described in detail 
mSi. The second group comprised 10 bicrystal SQUIDs, which were also single layer ones, using 150 - 200 nm thick 
YBaiCuzOi-x films deposited by linear hollow cathode sputtering onto SrTiOs bicrystal symmetrical substrates 
with 24° misorientation angleSI. All SQUIDs have a > 1 and T > 0.05 and all measurements were performed in liquid 
nitrogen at 77 K. 

Since the critical current, Ic of a high Tc DC SQUID cannot be measured with a good accuracy due to high level 
of thermal fluctuations, we here, for the comparison with experiment, use only the quantities, which arc measured 
directly, AVmax and I max- The dependence AVrmax {I max) for symmetric SQUID together with experimental 
points is shown on Fig. 1131 It can be seen, that there is a significant deviation of experimental points from theoretical 
line of symmetric DC SQUIDs. Part of the experimental points which lie above theoretical line can be explained by 
the junction asymmetry. In order to show this we add to the plot of Fig. ll3l the theoretical curves for asymmetric DC 
SQUID fFig. 114b and Fig. 114b). We choose a significant asymmetry on these graphs (7 = 0.9; p — 0.9) in order to 
mark the borders of the possible scattering of the voltage modulation values. In addition, the asymmetry allows one 
to explain the experimental points which lie below the symmetric line at small bias currents Imax < \QpA. Since 
at small bias currents Imax is close to Ic (see Eq. (|27|) ) this reduction of the voltage modulation is consistent with 
the result of section IIIC where we showed that at relatively small critical currents a junction asymmetry reduced the 
voltage modulation. Therefore, the junction asymmetry can explain the experimental values of the voltage modulation 
which lie above the line for symmetric DC SQUID and the points which lie below symmetric line at small bias currents. 
However, a significant reduction of the voltage modulation, which lies well below a symmetric line in Figs. Upl and 1141 
cannot be explained by the junction asymmetry. A possible explanation of these points is the presence of relatively 
large amplitude of a second harmonic in the junction current phase relation^. 

V. CONCLUSION 

In the paper we consistently compared the analytical theory for the voltage-current characteristics of the large 
inductance (L > lOOpH) high- Tc DC SQUIDs that has been developed previously^ with the computer simulations 
and the experiment. It is shown that the theoretical voltage modulation for symmetric junctions is in a good agreement 
with the results of known computer simulations. It is also shown that the asymmetry of the junctions results in the 
increase of the voltage modulation in the large range of critical currents and asymmetry parameters. We compared 
our theory with the experimental values of the voltage modulation. It appeared that the asymmetry can explain a 
large portion of experimental values of the voltage modulation which lie above the theoretical curve for symmetric 
DC SQUID. It also explains experimental points which lie below the curve at small critical currents. However, a 
significant portion of these values which lie below the curve cannot be explained by the junction asymmetry. From 
our opinion a possible explanation of these low lying points is the presence of relatively large amplitude of a second 
harmonic in the junction current phase relationiS. 
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FIG. 1: VCC of symmetric DC SQUID at zero magnetic field and different inductances. Ic = W/J.A, T = 77K, the curve with 
a = (black box) was calculated from expression Q, the curves with a — 1.5 (triangle) and a — 10 (cross) were calculated 
from the expression (1101 . 
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FIG. 2: The reduced voltage modulation vs bias current curves at different values of the critical current for symmetric DC 
SQUID. 
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FIG. 3: The maximum value of the reduced voltage modulation, AVr z max vs critical current. The different points on the 
curve correspond to the different bias currents I max, at which AV reaches its maximum. 
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FIG. 4: Dependence of the maximum value of the reduced voltage modulation AVr,max on the bias current Imax- 
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FIG. 6: The dependence of the SQUID critical current lo = 21c on the bias current I max- Solid line is the expression of 
Drung (Eq. ; crosses are theoretical points obtained from Eq. 1111 . 
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FIG. 7: The dependence of the normalized voltage modulation on the bias current. Comparison of theory (stars) with the 
simulation of Kleiner (opened circles). The curve for L = was calculated from J7J and JSJl, the curves for L — 94pT, L — 157pT 
were calculated from il l U . All calculations are made at 77 K for V = 1. 
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FIG. 9: The influence of the geometric asymmetry on the maximum voltage modulation, AVr,max- 
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FIG. 10: The dependence of the maximum voltage modulation, AVr,max on the critical current for two types of asymmetry: 
a) the current asymmetry (p = 0); b) the resistance asymmetry (7 = 0) 
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FIG. 11: The dependence of reduced maximum voltage modulation on the asymmetry parameters at Ic = 5pA for three types 
of asymmetry: a) p — 0, 7 =fc 0; b) 7 = 0, p 0; c) p — 7 7^ 0. 
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FIG. 12: The dependence of reduced maximum voltage modulation on the asymmetry parameters at Ic = 30pA for three types 
of asymmetry: a) p = 0, 7 / 0; b) 7 = 0, p 7^ 0; c) p = 7 0. 
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FIG. 13: The dependence AVr,max{Imax) for symmetric SQUID together with experimental points. 
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